Electronic thermal conductivity of disordered metals 
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We calculate the thermal conductivity of interacting electrons in disordered metals. In our anal- 
ysis we point out that the interaction affects thermal transport through two distinct mechanims, 
associated with quantum interference corrections and energy exchange of the quasi particles with 
the electromagnetic environment, respectively. The latter is seen to lead to a violation of the 
Wiedemann- Franz law. Our theory predicts a strong enhancement of the Lorenz ratio n/aT over 
the value which is predicted by the Wiedemann-Franz law, when the electrons encounter a large 
environmental impedance. 



PACS numbers: 72.15.Eb 

The Wiedemann-Franz law relates the electronic ther- 
mal conductivity, k, and the electrical conductivity, a, 
and states that the Lorenz ratio L = k/ctT is a univer- 
sal constant given by L = 7r 2 fcfj/3e 2 . In this equation 
ks is the Boltzmann constant, — e the electron charge, 
and T the temperature. The validity of the Wiedemann- 
Franz law relies mainly on a single-particle description of 
the transport properties, on the Fermi statistics of the 
charge carriers and on the assumption of purely elastic 
scattering 1 . In a Fermi- liquid, one expects that this law 
still holds at low enough temperature, when the quasi- 
particles cannot exchange energy during collisions 2 . De- 
viations from the Wiedemann-Franz law as recently ob- 
served in the normal state of a copper-oxide supercon- 
ductor have thus been interpreted as an evidence for the 
breakdown of Fermi-liquid theory^. 

The effects of Coulomb interaction on the electrical 
transport at low temperature can be, broadly, grouped in 
two main types. From one side, transport implies adding 
charges to a conductor. This has an energy cost that 
depends on the size and shape of the conductor itself. 
For example in tunnel junctions the energy transfer be- 
tween quasi-particles and the electrodynamical environ- 
ment causes the Coulomb blockade phenomena. On the 
other side Coulomb interaction leads in a disordered con- 
ductor to an additional source of random scattering that 
interferes with the scattering from the impurities. This 
is a quantum effect and depends on the details of charge 
diffusion. 

In the 1980s these quantum interferences were shown 
to lead to corrections to the electrical conductivity be- 
yond the standard Fermi-liquid results. It turned out 
that these corrections may, in fact, be incorporated 
into a scale-dependent renormalization of the Landau 
Fermi-liquid parametera^SiLSiSi. Whithin this frame- 
work, Castellani and co-workers^ demonstrated that the 
scale-dependent corrections to the thermal conductivity 
are the same as the corrections to the electrical conduc- 
tivity. This led them to conclude that the Wiedemann- 
Franz law is valid up to the metal-insulator transition. 



This conclusion was challenged by Livanov et al^i 
who found for a two-dimensional system with long range 
Coulomb interaction additional contributions to the ther- 
mal conductivity. Recently, the issue has been reexam- 
ined by Niven and Smith 12 , who also concluded that the 
Wiedemann-Franz law is violated. 

In this paper we study the problem by means of 
the quasiclassical Green's function approach, which has 
proved to be a powerful tool in describing the dynam- 
ical properties of superconductors^ and the transport 
in hybrid mesoscopic structures^. Recently it was fur- 
ther demonstrated that both Coulomb blockade phenom- 
ena and quantum interference corrections to the charge 
transport can be conveniently described within this theo- 
retical framework^. Advantages of the method are that 
it is not restricted to the linear response regime, and of- 
ten provides more compact derivations than the standard 
diagrammatic techniques. 

A perturbative calculation of the thermal conductivity, 
besides confirming Ref. 12 , allows us to clarify the origin 
of the apparent discrepancies in the literature. To do so 
we separate the different physical mechanisms by their 
different range of exchanged energies and relevant length 
scales. For instance, the quantum interference effects oc- 
cur over distances from the mean free path up to the 
thermal diffusion length \J HD/ksT, and imply energy 
exchanges larger than the temperature T. These yield 
corrections which are logarithmically divergent and can 
be readily related to the scale-dependent renormalization 
of the electrical conductivity. Here the temperature acts 
as an infrared cutoff. The interaction effects responsi- 
ble for the deviations from the Wiedemann-Franz law 
are associated with the long-range part of the Coulomb 
interaction and their singular behavior has the tempera- 
ture as the upper cutoff. We successively concentrate on 
the long-range part of the Coulomb interaction. In par- 
ticular we predict a sizable enhancement of the Lorenz 
ratio when the sheet resistance is of the order h/e 2 in the 
case of a two-dimensional electron system, and a strong 
enhancement of the Lorenz ratio for thin metallic wires 
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when the total resistance of the wire is larger than h/e 2 . 

From now on we use units where h = kg = 1, but we 
put back the constants in the final results. We start with 
a brief introduction of the quasiclassical formalism. For 
a more detailed description we defer the reader to Rcf. 16 . 
Whithin this formalism the short-distance behavior of 
the electron Green's function is taken into account in an 
averaged way from the outset by introducing the quasi- 
classical Green's function, which solves the Eilenberger 
equation^ 

[d tl + d t2 + v F p ■ <9 X ] g tl t 2 (x, p) = -i [S (x, p) , g (x, p)] . 

In contrast to the Dyson equation for the ordinary 
Green's function, the Eilenberger equation for g is homo- 
geneous and requires a normalization condition, which 
can be chosen of the form gg = 1. The Green's function 
has a two-by-two matrix structure in Keldysh space, 
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g R g K 

g A 



(2) 



Matrix products imply both summation and integration 
over Keldysh indices and time variables, respectively. 
We recall that, whereas the diagonal components of g 
describe the spectral properties of the system, the off- 
diagonal Keldysh component carries information about 
the distribution function. In this respect, the Keldysh 
component of Eq. Q is the quantum analog of the Boltz- 
mann equation. 

Impurity scattering is introduced by means of the stan- 
dard white-noise random potential and is described by 
the self-energy in the self-consistent Born-approximation 
as 

i f dp , „ , 
s tit 2 ( x ) = -^: / ^<?ti*2 ( x - p) , (3) 

where r is the elastic scattering time and fid is the d- 
dimensional solid angle. 

The charge and heat current densities have the form 



j Q (x,i) / 2 * J fid F ' 



9 (x, t; p, e) , 
(4) 

where No is the density of states at the Fermi energy, 
t = (ti+t2)/2 and e corresponds to the Fourier transform 
of the relative time t\ — t2- 

In the dirty limit, the variation of the Green's function 
is on space and time scales larger than the elastic mean 
free path I = vpT and scattering time r, respectively. 
In this limit one may expand the Green's function g in 
spherical harmonics and keep the s- and p-wave compo- 
nents only, g(p) — g s + pg p + ... The Eilenberger equation 
is then replaced by ( D = vpr/d) the Usadel equation^, 



d t g s {x-) - Dd x [g s d x g s ] = 



(5) 



which is the analogous of the saddle-point condition in 
the non-linear cr-model matrix field theory^. As a re- 
sult, the currents are expressed in terms of the s-wave 



component of the Green's function: 

j e (x,t) = - eN ° D J de(g s d x g s ) K , (6) 

JqCM) = ^tt- [ dee (g s d x g s ) K . (7) 



As a simple application of the formalism we derive the 
Drude formula for the electrical and thermal conductiv- 
ity. In the absence of interactions the Green's function g 
reads: 



§ s (x,t;e) 



1 2F 
-1 



(8) 



Near local equilibrium with a local temperature T(x) and 
chemical potential /i(x) the function F is given by 



F = tanh 



£ - M(x) 
2T (x) J ' 



(9) 



from which the Drude expressions for both electrical and 
heat currents are found 



j e = 2e 2 DN Q (Vfi/e) 
j Q = ^-k%2N DT(-VT), 



(10) 
(11) 



and in particular the Wiedemann-Franz law holds. 

To include the effects of Coulomb interaction, we 
introduce 15 a Hubbard-Stratonovich matrix field 



i>2 4>x 



(12) 



whose fluctuations describe the retarded, advanced, and 
Keldysh components of the screened Coulomb interaction 
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In the presence of the field </>, one first adds a term ie[</>, g] 
to the right-hand side of the Eilenberger Q or the Us- 
adel equation © ■ Secondly, the resulting solution g [0] is 
averaged over the fluctuations of 4> according to Eq. i|13[l • 
In analogy to the non- interacting case (Cf. Eq. JSJ)), one 
can define the distribution function in the presence of in- 
teractions via the relation between the Keldysh and the 
retarded, advanced components of the Green's function 
(g K ) = (g R )F — F(g A ). We further assume a system 
which is - with the exception of a weak temperature gra- 
dient - translational invariant. Then it is convenient to 
expand the distribution function as 

F £ _ w ( Xl ) « F £ _ w (x) - d T F e _ w (x)VT • ( Xl - x) (14) 

and to Fourier transform from real to momentum space. 
The correction to the thermal current is finally obtained 
as Sj Q = <5j% + Sj b Q with 

6j a Q = DN Q VT J dee J ^9 T (^-c( x )F e (x)) 

rK?(q) (15) 
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and 



DN VT J dee J ^F e (x)d T F £ _ w (x) 
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-iw + Dq 



2Y>> 



Kf(q), (is) 



where d is the dimension of the system under consider- 
ation. Our result, Eqs. 1(15116(1 . for the thermal current 
is equivalent to the thermal conductivity found in»£ by 
using the diagrammatic method and the Matsubara tech- 
nique. We notice that the diffusive pole appearing in 
Eqs. (|f 5116(1 originates from the Usadel equation JSJl. 

Using the relation F e F £ _ w = 1 - (F e - F e ^)B{uj /2T) 
with B(x) = coth(x) allows to evaluate the e-integrations 
in Eqs. 1(15116(1 with the result 



deed T (F e F e _ u ) = -oj 2 d T B (^) , 
2vr 2 T 

deeF e d T F,- u = —d^uB [ , 



(17) 



3T 0ul \2T 



(18) 



From Eqs. (|17|1 and l|18() we observe that 5jq is dominated 
by diffusive modes of frequency \u>\ < T, whereas modes 
with frequencies |u>| > T give the dominant contribution 
to <5jg. To appreciate the role played by the different fre- 
quency ranges we begin by evaluating the current in two 
dimensions. The retarded component of the dynamically 
screened Coulomb interaction reads 



V R (q, W ) 



1 K 2 d — iw + Dq 2 
2N Q q -ioj + Dn 2 dq ' 



(19) 



where K2d = 47re 2 ./Vo is the screening vector in two di- 
mensions. By considering first <5jg one notices that the 
momentum integration to be performed is identical to 
the momentum integral in the correction to the electrical 
conductivity^^, i.e., the integration is logarithmically 
divergent in the ultraviolet and must be cutoff with the 
diffusive condition Dq 2 r < 1. In the w-integration there 
is a minor difference at low frequencies \u>\ < T due to the 
second term on the right-hand side of Eq. 1(18(1 . In two 
dimensions, with logarithmic accuracy, this difference is 
negligible and one has 



t? M-vt), 



(20) 



where 5a = — e 2 /(2-K 2 ) ln(l/TY) is the interaction cor- 
rection to the electrical conductivity and Tt < 1. The 
other contribution to the thermal current, <5jg, does not 
depend on the ultraviolet cutoff 1/r, 
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since the temperature acts as an upper cutoff in the fre- 
quency integration. In contrast, in the limit of good 
metallic screening when K2d — * oo, the integration be- 
comes infared divergent. By combining the two contri- 
butions we finally write the expression for the thermal 
conductivity in a form which shows that, although the 
Wiedemann-Franz law is violated, 
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/k B T)\, (22) 



the integration of diffusive modes in the region T < 



Dq 2 ,ui 
and k, 



< t yields the same scaling equations for a 



din a d In k 



din I din I ' 



(23) 



so that the apparent discrepancies in the literature are 
no contradiction. 

We observe that in the last term of Eq. (|22f) . respon- 
sible for the violation of the Wiedemann-Franz law, only 
the extreme long wavelength modes of the dynamically 
screened Coulomb interaction with Dq 2 < \u\ < T are 
relevant, cf. Eqs. 1(1 5(1 and 1(17(1 . It has been shown in 
Refill that these can be summed to all orders and in the 
end modify the Green's function like a gauge factor, 



<?*it 2 ( x ;p) 



e i ^ x ' tl) P t lt2 (x;p)e- 



■il,5(x,t2) 



(24) 



where (9t</?(x, t) — e^(x, £). Whereas the gauge factors 
drop in the expression for the electrical current, i.e. the 
long wavelength modes of the Coulomb interaction do 
not modify the electrical conductivity, they survive in 
the heat current to yield 

JQ = J f^ VF P J de { e 9o (x,*;p,e) 

- \z ($(t)$(x, t; P, e) + 5(x, t; p, e)0(t)» K }. (25) 

Eq. 1(25(1 makes clear the physical origin of the violation 
of the Wiedemann-Franz law. While the first term on 
the right-hand side reproduces the non-interacting con- 
tribution to the thermal current, the second may be in- 
terpreted as the effect of the time dependent fluctuations 
of the quasi-particle energy in the presence of an elec- 
tromagnetic environment. Indeed, the extra heat cur- 
rent is proportional to the correlation of voltage and cur- 
rent fluctuations in the system, 5]q — (</>i(x, t)<5j e (x, t)), 
which then leads to the strikingly simple result 

JQ = y^f a(-VT) - ^V(M^t)M^t)). (26) 

Notice that due to the linear current- voltage characteris- 
tics of the system under consideration only the first order 
in the Coulomb interaction contributes to the heat cur- 
rent. By using the fluctuation dissipation theorem (or 



4 



equivalently Eq. (JT3J)) 



<0i(x,t)0i(x,t)) = 




direct contact can be made with the conventional pertur- 
bation theory, i.e. with 5]q in Eq. JTSJl. 

Instead of parameterizing the local voltage fluctuations 
in terms of an interaction V R (q, ui) we will in the follow- 
ing parameterize them in terms of the impedance of the 
local electromagnetic environment, Z(u>). By doing so, 
the thermal conductivity reads 

_n 2 k 2 B T ak B f f E/2k B T \ 2 RcZ(E/h) 

where for clarity we put back the H and k B . We will 
now discuss three different examples for the impedance 
Z . The simplest situation consists of a purely ohmic en- 
vironment where Z(E/h) = R. The thermal conductiv- 
ity is found linear in the temperature, strong deviations 
from the Wiedemann Franz law are found when the en- 
vironmental resistance is of the order of the resistance 
quantum h/e 2 or larger. The explicit result is 

3 e z 

From the retarded Coulomb interaction as given in 
Eq. (|19J) we determine the impedance of a thin film as 
ReZ{E/h) = {l/4ncr)\n(hDK 2 d k B T/E 2 ). Due to the 
weak logarithmic energy dependence of the impedance 
the thermal conductivity is to good accuracy obtained 
from Eq. © with R = ReZ(k B T/h). As a third ex- 
ample we consider a i?C-transmission line, as a model of 
a gated wire. The impedance is TLeZ = ^y/Ro/2\u>\Co, 



where Rq and Co are the resistance and capacitance per 
unit length. We find a contribution to the thermal con- 
ductivity which is proportional to the square root of the 
temperature, 



K = ~~3~e 2 ~ <J 2 - 456 ^ ff hp ' (30) 

with 2.456 the approximate numerical value for 
3C(3/2)r(3/2)/2 3 / 2 . 

In summary we calculated the thermal conductivity 
of disordered metals. In the two-dimensional electron 
system the scaling equations for the thermal and the 
electrical conductivity are the same, nevertheless the 
Wiedemann-Franz law does not hold. The deviations 
from the Wiedemann-Franz law are comparable in size 
to the localization effects. It is interesting to note that 
this is in qualitative agreement with observations made in 
the cuprate o 3 ' 20 : The resistivity of PCCO in Refi 3 - shows 
a well pronounced low temperature anomaly which has 
been attributed to localization effects, and at the same 
time the low temperature heat conductivity is larger 
than what would be expected from the Wiedemann- 
Franz law. In the low temperature resistivity of Tl- 
2201 in Ref. 20 no indications of localization effects are 
seen and the Wiedemann-Franz law is perfectly obeyed 
within the experimental accuracy. Quantitatively on the 
other hand the agreement of our theory with Ref^ re- 
mains poor, since the sheet resistance was estimated as 
ps h/(60e 2 ) from which we expect a much smaller 
enhancement of the heat conductivity than observed ex- 
perimentally. By measuring the Lorenz ratio in a gated 
film or wire as function of the gate capacitance, it should 
be possible to test our predictions experimentally. 
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